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ABSTRACT 

Instabilities in planetary systems can result in the ejection of planets from their host 
system, resulting in free-floating planets (FFPs). If this occurs in a star cluster, the 
FFP may remain bound to the star cluster for some time and interact with the other 
cluster members until it is ejected. Here, we use IV-body simulations to characterise 
close star-planet and planet-planet encounters and the dynamical fate of the FFP pop¬ 
ulation in star clusters containing 500 — 2000 single or binary star members. We find 
that FFPs ejected from their planetary system at low velocities typically leave the star 
cluster 40% earlier than their host stars, and experience tens of close (< 1000 AU) 
encounters with other stars and planets before they escape. The fraction of FFPs 
that experiences a close encounter depends on both the stellar density and the initial 
velocity distribution of the FFPs. Approximately half of the close encounters occur 
within the first 30 Myr, and only 10% occur after 100 Myr. The periastron veloc¬ 
ity distribution for all encounters is well-described by a modified Maxwell-Bolzmann 
distribution, and the periastron distance distribution is linear over almost the entire 
range of distances considered, and flattens off for very close encounters due to strong 
gravitational focusing. Close encounters with FFPs can perturb existing planetary 
systems and their debris structures, and they can result in re-capture of FFPs. In 
addition, these FFP populations may be observed in young star clusters in imaging 
surveys; a comparison between observations and dynamical predictions may provide 
clues to the early phases of stellar and planetary dynamics in star clusters. 

Key words: Open clusters and associations: general; stars: kinematics and dynamics; 
planets: dynamical evolution and stability 


1 INTRODUCTION 


The majority of stars form in clustered environments (e.g., 
Lada & Lada 2003 Portegies Zwart et al.|20l0 l, where close 
encounters between stars are frequent. These close encoun¬ 


ters can perturb or destroy planetary systems (e.g., Spurzem 
et al.|2009 Boley et al.|2012 I. Even a mild perturbation can 
break up a marginally stable planetary system and can even 


affect the shortest-period planets in the system (e.g., Hao et 
which can result in strong planet-planet scatter- 


al. 


2013) 


ing, physical collisions between planets and their host stars 


2008 Juric & Tremaine 

2008 Nagasawa & Ida]2011 

1, long- 

term secular evolution (e.g., Malmberg et al. 2007 

Parker 




* E-mail: long.wang@pku.edu.cn 


et al. |2009 Malmberg et al.|2011 1, and the ejection of plan¬ 


ets from the system. In addition, mass loss due to stellar 
evolution of single or binary host stars can result in loss 
of planetary companions through similar interactions (e.g., 


Veras et al. 2011. Veras & Tout||2012; Adams et al. 2013; 

Voyatzis et al. 

2013 

Nowak et al. 

2013 

). Even the Galac- 


tic tidal field can indirectly play a role in the disruption of 
planetary systems when a wide stellar companion is present 
(e.g., Kaib et al. 20131. Veras & Raymond (20121 demon¬ 


strate that the dynamical and stellar evolution of isolated 
planetary systems alone cannot account for the observed 
free-floating planet (FFP) population in the Galactic field, 
and that close encounters in star clusters are an important 
source of FFPs in the field. 

Although only a small number of FFP candidates have 
been detected so far, they are potentially abundant in the 
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Galactic disk (e.g., Strigari et al. 2012). Exoplanets that Table 1. Initial conditions for the modelled star clusters. 


are too distant from stars to directly affect their observable 
properties are most easily detected using microlensing (e.g., 
Mao fe Paczynski|1991||Gould fe Loeb|1992||Abe et al.|2004} 


Quantity 


Value 


Beaulieu et al.|2006 Gaudi|2012 1. Microlensing surveys also 

have the potential to discover FFPs (e.g., Di Stefano|2012 1, 
and results imply that there are roughly twice as many FFPs 
in the Solar neighbourhood than there are main-sequence 
stars (Sumi et al7|[2C)l 1 ). In addition, deep imaging surveys 
can be used to detect young planetary-mass objects near the 


Number of stars & binaries 
Half-mass radius 
Dynamical model 
Virial ratio 
Tidal field 


N = N s + N b = 500, 1000, 2000 
r hm = 0.38, 0.77, 1.54 pc 


(1911) 


Delorme et al. 


et al. 2006 

Caballero et al. 2007; Bihain et al. 2009; Quanz 

et al. 2010 

Pena Ramirez et al.| 2011; jScholz et al. 2012 


20121. These observational studies help us 


further constrain the origin and fate of these FFPs provided 
that a good understanding of their dynamics is known. 

FFPs are thought to have been ejected from their sys¬ 
tem with velocities of typically 0.1 — 10 kms -1 , as a result 
of planet-planet scattering (delayed ejection) or immediately 
after a close encounter with a passing star (prompt ejection), 
and the ejection velocities from the latter process tend to be 
higher (Malmberg et al. 2011| ). When a FFP is ejected in 
a star cluster, it may escape immediately if its ejection ve¬ 
locity exceeds the local escape velocity of the star cluster, 
or it may remain bound to its host cluster for millions of 
years, until it escapes through ejection or evaporation. Dur¬ 
ing its life in a star cluster, a FFP can experience multiple 
close encounters with other stars before escaping, and may 


Initial mass function 
Binary fraction 
Semi-major axis distr. 
Mass ratio distr. 
Eccentricity distr. 
Orbital orientation 


Plummer 

= 1/2 

Galactic Solar orbit 

Kroupa et al. (1993 I, 0.2 — 5 Mq 
T3 = NjJN = 0%T20%, 50% 
faip.) oc a -1 ; (10 -6 - 10 -3 ) X r vi 
Random pairing 
/ e (e) = 2e (0 < e < 1) 

Random 


Planet-to-star ratio 
Planet mass 
Density distribution 
Planet ejection velocity 


77 = N p /N sh = 0.5,1,2 
M p = 1 Mj 

Equivalent to that of stars 
M s , f e j (see Figure |T] and § 2.2 i 


and lV a b = N a + 2N b is the total number of stars in the clus¬ 
ter. The initial values for the virial radii are r v i r = 0.5, 1.0, 
and 2.0 pc, which are typical for open clusters (e.g. 


Lada 


& Lad a|2003 1, and correspond to intrinsic half-mass radii of 
'rhm = 0.38, 0.77 and 1.54 pc, since for the Plummer model 


1.30r hm (e.g., Heggie & Hut 2003). Each star clus 


even be re-captured by another star l 

Kouwenhoven et al. 

2010; Malmberg et al. 2011. Moeckel & Clarke 2011, Parker 

& Quanz 

2012 

Perets & Kouwenhoven| 

2012). Direct N- 


body simulations of single-planet systems in star clusters 
have shown that many are disrupted and that the resulting 
FFPs can remain in these star clusters for many millions 
of years (e.g., Hurley fe Shara||2002 Parker fe Quanz]|2012 


Craig & Krumholz 20131. Fly-by simulations mimicking the 


evolution of multi-planet systems in star clusters confirm 
that the survival of these systems depends strongly on the 
properties of the stellar environment and the semi-major 
axes of the planets, but also demonstrated that planetary 


ter is initially in virial equilibrium, i.e., Q = \K/P\ = 1/2, 
where K and P are the total kinetic and potential ener¬ 
gies, respectively. As a consequence, the corresponding ini¬ 
tial (one-dimensional) velocity dispersion a(r) as a function 
of distance to the cluster centre r is 

2 , , GM / r 2 \“ 1/2 

where G the gravitational constant and M the total cluster 
mass ( Heggie fe Hut||2003 l. 


Following Malmberg et al. 


drawn from the Kroupa et al. 
(IMF) in the mass range 0.2 — 5 


(2007), the stellar masses are 


(1993) initial mass function 
There is evidence 


Afm. 


turbed systems also mutually interact (e.g., Chatterjee et al. 

2012; Hao et al. 2013; Liu et al. 2013). 

star cluster (e.g., Weidner & Kroupa 2004; Weidner et al. 

2013), this dependence is still under debate (e.g., Cervino 

The aim of this study is to analyse the dynamical prop- 

et al.|2013a|b). We therefore follow 

Malmberg et al. (2007) 


focus on close encounters between the members of the star 
cluster (single stars, binary stars, and FFPs). This article 
is organised as follows. The methods and assumptions are 
described in §[2] The results are presented in §|3] Finally, we 
draw the conclusions and describe our future work in § [4] 


2 METHOD AND ASSUMPTIONS 

2.1 Initial conditions for the star clusters 

The properties of our model star clusters are summarised in 
Table [l] We create star clusters with [Plummer (1911) den¬ 
sity and velocity distributions. We carry out simulations of 
open star clusters with N = N s + N b = 500, 1000 and 2000 
systems, where N s and N b represent the number of single 
stars and binary systems in the star cluster, respectively, 


in choosing a constant upper mass limit, while realising the 
possibility that realistic clusters may host a more massive 
star that can alter the dynamics of the FFP population. 

We carry out the simulations with binary fractions B = 
N b /N = N h /(Ns+N h ) = 0%, 20%, and 50%. The individual 
components of the binary systems are randomly paired from 
the IMF (see Kouwenho ven[ 2006 Kouwenhoven et al. 2009, 
for details). For this method of pairing stars into binary 
systems the total mass of the star cluster equals N(M)( 1 + 
B), where (M) is the average stellar mass. The adopted semi¬ 
major axis distribution is f a (a ) oc a -1 , which corresponds 
to a flat distribution in logo, also known as Opik’s law (e.g., 


van Albada|1968. Vereshchagin 

eTaLll987; 

Poveda & Allen 

2004j Kouwenhoven et al.|2005[ 

2007). Sem 

i-major axes are 
“ 6 —10 -3 ) xr vir , 

drawn from this distribution in t 

re range (1C 


which corresponds to 1 — 1000 AU for the star clusters with 
rhm = 0.38 pc. We do not include very tight binary systems 
as they are relatively inert and effectively act as single stars. 
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Figure 1. The cumulative distribution of velocities at which plan¬ 
ets escape from their host stars, as obtained by Malmbcrg ct al. 
(2011). The four curves indicate the results for models with differ¬ 
ent host star masses, and represent the distributions for prompt 
ejections (during the encounter) and delayed ejections (after the 
encounter). For comparison, we also model the evolution of FFP 
populations with zero ejection velocities, in addition to the four 
velocity distributions shown here. 


We also do not include very wide binary systems, as their 
individual components experience encounters very similar to 
those of single stars. Moreover, many of these wide binary 


systems are easily destroyed (e.g., Parker et al. 20091. We 
choose a thermal eccentricity distribution f e (e) = 2e for 0 < 
e < 1 (Heggie 1975), and all orbits are assigned a random 


spatial orientation and a random orbital phase. 

We include an external tidal field by modelling the star 
clusters on a Galactic circular orbit in the solar neighbour¬ 
hood. The tidal radii of the star clusters are then given by 
1/3 / \ 1/3 


n(M) = 


GM 


4A(A - B) 


1.41 


M 

Mo 


pc , 


( 2 ) 


where A = 14.4kms _1 kpc _1 and B = — 12.0kms _1 kpc _1 
are the Oort constants (e.g., Binney fe Tremaine|1987 l, and 
M is the total mass of the cluster. This external tidal field 
enhances the dissolution of the cluster, in particular for low- 
mass and low-density star clusters. In our grid of models the 
tidal radii are rt = 6.0 — 11.5 pc, depending on total mass 
of the star cluster. 


2.2 Initial conditions for the free-floating planets 

In addition to the stellar population we also add a popula¬ 
tion of FFPs. The initial positions of these planets are drawn 
from the same distribution as that of the stars. Since FFPs 
are thought to have been ejected from planetary systems, 
they have an initial position distribution which is statisti¬ 
cally identical to that of their host stars. It can be argued 
that planets are preferably ejected from their host systems 
near the cluster centre, as the close encounter frequency with 


other stars is significantly larger in this region. On the other 
hand, from delayed ejection in perturbed multi-planet sys¬ 
tems, we can expect an initial spatial distribution of ejected 
planets that is similar to that of the stars in the star clus¬ 
ter when the delay time is larger than the crossing time, 
which is 1 — 3 Myr for our cluster models. Fly-by experi¬ 


ments by Malmberg et al. (20111, Hao et al. (20131 and Liu 


et al. (2013) have demonstrated that, although some close 


encounters result in prompt ejections, perturbed planetary 
systems can eject planets up to tens of millions of years af¬ 
ter a close encounter has taken place, at a location that is 
unrelated to the place where the close encounter took place. 

Moreover, we make the assumption that all FFPs are 
free-floating at the start of our simulations, and that none 
of the stars have planetary companions. Hence, we ignore 
the possibility of new FFPs being ejected from perturbed 
planetary systems in the star clusters. Although this may 
seem a strong constraint to the applicability of our study, we 
show in Section[3]that many of our results are also applicable 
to populations of FFPs that are ejected at different times. 
Moreover, our results can also be use to calculate encounter 
probabilities for individual FFPs, and to estimate the time 
at which an individual FFPs may escape. 

Two dynamical limits can be considered here: (i) the 
limit where all planets have zero ejection velocity with re¬ 
spect to their host stars, and (ii) the limit where all planets 
have ejection velocities that are substantially larger than 
the star cluster’s escape velocity. The latter case is triv¬ 
ial and results in the immediate ejection of all FFPs. In 
the former case, both the stellar and planetary populations 
initially have identical position and velocity distributions. 
Subsequently, energy exchange between the bodies, result 
in mass segregation and also in the gradual ejection of all 
planets from the system. 

The adopted ejection velocity distributions f e j are based 
on the work of |Malmberg et al.| ( |2011| ), who studied the 
evolution of planetary systems consisting of four gas-giant 
planets orbiting a star in a star cluster. They studied the 
prompt ejection velocities of planets escaping during close 
encounters with other stars (/ e j = P), and the delayed ejec¬ 
tions as a result of planet-planet scattering in a perturbed 
system (/ e j = D). Their simulations were carried out using 
host stars with masses of M s = 0.6 Mq and M s = I.OMq. 
The cumulative ejection velocities for the different models 
are shown in Figure [T] For comparison, we also carry out 
simulations with zero ejection velocities, f e j = Z, where the 
initial velocity distributions of planets and stars are iden¬ 
tical. We study the dynamical evolution of the FFP pop¬ 
ulation for each of these five distributions / e j. As the two 
delayed ejection velocity distributions in Figure [T] are simi¬ 
lar, the results for the models with M s = I.OMq and f e j = D 
are omitted in most figures and discussion throughout the 
remaining part of this paper. 

Each planet is assigned a mass of Mp = 1 Mj as 
1O -3 M0. Since these planets are of very low mass compared 
to the stars, the actual value of the mass does not affect 
the evolution of the stellar population and star cluster as a 
whole, the star-planet encounter rate, or the rate of planet 
capture by stars. It may, however, slightly affect the planet- 
planet encounter properties. 

For the same reason, the total number of FFPs does 
not affect the evolution of the stellar population in the star 
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cluster. It simply scales the total number of close encounters 
and capture events involving planets. Inspired by the find¬ 


ings of Sumi et al. (2011) we adopt an initial planet-to-star 
ratio of 77(0) = N p (0) /IV a b(0) = 2 in our default model, and 
also carry out simulations with 77(0) = 0.5, and 77(0) = 1, 
for comparison. 

As planets are presumed to have been ejected from the 
planetary system in which they formed, the velocity disper¬ 
sion of the FFP population is larger than that of the stel¬ 
lar population. The initial velocity of each FFP therefore 
constructed by combining two velocity vectors. The main 
component, v P (r), representing the velocity of the host star 
at the moment of ejection, is drawn from the the Plummer 
distribution (see Eq.[l]). The second velocity component, v e , 
is drawn from the distribution of ejection velocities at which 
the FFP escaped from its host star: 


UFFp(r) = Vp + Ve ■ 


( 3 ) 


The resulting initial velocity dispersion for the FFP popu¬ 
lation is then 

ctffp = Op(r) + a\ , (4) 

because we draw v p and v e independently and with random 
spatial orientations. 

The (one-dimensional) velocity dispersion of the stars 
at the half-mass radius (Eq. [l]) ranges between fr p (rh m ) « 
0.3 kms" (for 7V ab = 500 and r hm = 1.54 pc) and cr p (rh m ) ~ 
1.5 kms -1 (for N sb = 3000 and rhm = 0.38 pc). The planet 
ejection velocities are typically in the range 0.1 — 5 kms -1 , 
and can therefore contribute substantially to the velocity 
dispersion of the FFPs in the star clusters. We therefore 
expect a relatively large fraction of the FFPs to escape at 
early times from low-density star clusters. 


2.3 Simulation parameters 


We use the NBODY6 package ( |Aarseth| 1 999| 2003) to carry 
out the simulations. Stellar evolution is modelled following 


the prescriptions of Eggleton et al. (1989 19901 and Hurley 
|et al. | (|2000|). The NBODY6 code was modified to suit the 


requirements of (i) the inclusion of the free-floating planets, 
and (ii) identification of the close encounters. 

In a star cluster with a virial radius r vlT consisting of 
N single stars, when a pair of stars approaches within a 
distance r c i oae « 4r v i T /N their velocities are significantly 
perturbed (e.g., |Aarseth] [2003)>. In our grid of simulations, 
J'dose is in most cases several hundred to somewhat more 
than a thousand astronomical units. Inspired by the previ¬ 
ous work of |Malmberg et al.~| ( |2007| ) , we record encounters 
within a distance of r en c < 1000 AU. Encounters with hard 
(i.e., tight) binaries are registered as such, while for encoun¬ 
ters with soft (i.e., wide) binary systems, the encounters 
with individual stars are registered as encounters with sin¬ 
gle stars. The individual components of binary systems with 
periastron separations smaller than 1000 AU experience en¬ 
counters on a regular basis, according to this definition, but 
are excluded from the analysis of the close encounter prop¬ 
erties. All simulations are carried out until the star cluster 
has completely dissolved, such that we are able to register 
all close encounters that occur during the lifetime of each 
star cluster. In order to account for statistical effects, we 
simulate ten realisations for each model. 




Figure 2. Top: the evolution of the number of cluster member 
stars N sb (blue curve) and planets A), (red curve) as a function of 
time. This particular example shows the results for the model with 
N s b = 2000, B = 0%, 77(0) = 1, and rb m = 0.38 pc. The distri¬ 
bution at which the FFPs are assumed to have been ejected from 
their host star is identical to the model with f e j = P (for 1 A7q, ) 
in Figure]!] Bottom: the planet-to-star ratio 77(f) = N p {t)/N sb (t) 
as a function of time t. The red line indicates the best linear fit 
(Eq.[l5). 


3 RESULTS 


3.1 Star cluster evolution and membership 

3.1.1 Star cluster dynamics 

The dynamical evolution of the stellar population in the star 
cluster is practically unaffected by the presence of the plan¬ 
ets, as the planets are three orders of magnitude lower in 
mass than the stars, and therefore effectively behave as test 
particles. The evolution of the stellar and planetary popula¬ 
tions over time is illustrated in Figure [2] which shows that 
the number of FFP members in the cluster decreases faster 
than the number of stellar members. The reason for this 
twofold: (i) the FFP population is initially to some degree 
supervirial (see Eq. [4]), which results in the rapid escape 
of many of the FFPs with a high initial velocity, and (ii) 
dynamical interactions result in FFPs obtaining higher ve¬ 
locities than stellar components, resulting in the preferential 
loss of planets. 

The initial escape velocity v esc for the Plummer model 
as a function of distance to the cluster centre r is 


=(r) = ( 


2 GM\ 1/2 


1 + 


r 2 

r hm 


- 3/4 


( 5 ) 
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Figure 3. The fitted linear parameters IZ f (0) vs. af for planet- 
to-star ratio as a function of time (Eq. |15| ) for all models. Top : 
fixed r^m an d ejection velocity distribution of planets, for clusters 
with different AT s h = N s + and N p . Middle: fixed N p and ejec¬ 
tion velocity distribution of planets, for different rh m and N^. 
Bottom : fixed rh m and N p , for different planet ejection velocity 
distributions and different AT s h. The adopted initial planet ejec¬ 
tion velocity distributions include prompt ejection (P), delayed 
ejection ( D ), and zero ejection velocity (Z); see Figure [l] The 
symbols and their colours represent the different initial conditions 
for the simulations, as indicated in the legends. Each symbol rep¬ 
resents the average result for an ensemble of star clusters. 



Figure 4. The relation between the half-number time thp °f the 
number of planets in the clusters and the half-number time £hs 
of the number of stars in the clusters, for different sets of star 
clusters models and ejection velocity distributions. All models 
have N p = 1000. The three panels show the results for half-mass 
radii of rh m = 0.38, 0.77, and 1.54 pc. Each symbol represents the 
average result for an ensemble of ten star clusters with a certain 
N s b = N s ~h Nb- The colours indicate the different initial planet 
ejection velocity distributions for prompt ejection (P), delayed 
ejection ( D ), and zero ejection velocity (Z); see also Figure [l] 
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(Heggie & H ut|2 0031, where M is the total cluster mass and 
G the gravitational constant. At the half-mass radius (r = 
r bm ) this equation reduces to v eS c ~ 0.84-^/GM/r bm . For 
our models its value typically ranges between v eS c(rhm) ~ 
0.7 kms^ 1 (N a b = 500 and rh m = 1.54 pc) and rw^hm) ~ 
3.6 kms -1 (7V sb = 3000 and r bm = 0.38 pc). Although a 
number of FFPs are able to escape the cluster immediately, 
the FFPs with a relatively small ratio a e /a p (r) are able 
to remain part of the star cluster for some time, but are 
typically removed at earlier times than the stars. The stel¬ 
lar population is barely affected by these encounters, but 
the typical velocity of the FFPs in the star cluster rapidly 
increases. Two-body encounters typically occur on a relax¬ 
ation time. For stars with a mass close to the average stellar 
mass in a Plummer sphere, the (initial) relaxation time f bm 
at the half-mass radius is 


thm — 


0.206iV sb r^ 2 


( 6 ) 


where N s b is the number of stellar-mass objects, M = 
N s h{M s )(l + B) the total cluster mass, and InA « lniV sb 
is the Coulomb logarithm (Binney & Tremaine 1987; Hcg- 
gie fc Hut||2003 l. The half-mass relaxation time in our sim¬ 
ulations ranges from f bm ~ 7 Myr (for 7V sb = 3000 and 
rhm = 0.38 pc) to thm ~ 30 Myr (for 7V sb = 500 and 
rhm = 1.54 pc). 

Encounters between the star cluster members also result 
in energy equipartition, also roughly on a timescale 


t ms — 


(M sb ) 

M max 


thn 


(7) 


(e.g., Spitzer 19871 where (M sh ) = (M s }( 1 + B ) is the aver¬ 
age mass of a particle participating in the mass segregation 
process (a single star or a binary star) and M max is the mass 
of the most massive particle in the cluster. In our simula¬ 
tions, (M s ) /M max ~ 0.1, and the mass segregation timescale 
therefore ranges from 0.3 Myr to 3 Myr for the set of models 
studied in this paper (but may be shorter if the star clus¬ 
ters are initially substructured, see Allison et al. 2009a|b I. 
These timescales, as well as the dissolution time of the star 
clusters, also depend on the range of stellar masses present 
in the cluster (e.g., Kouwenhoven et al.|2014 |. FFPs do not 
contribute to the mass segregation process in the sense that 
they do not affect the stellar population. However, they do 
obtain higher velocities as a result of close encounters with 
stars, also roughly on a timescale thm- 


3.1.2 Analytic estimates for the survival times 

The FFP population can be separated into three categories: 
(i) FFPs that escape the cluster immediately, (ii) FFPs that 
are ejected from the cluster as a result of strong three-body 
encounters, and (iii) FFPs that evaporate from the cluster 
following a series of weak encounters and an interaction with 
the Galactic tidal field. Whether or not a FFP can escape 
from the star cluster immediately depends on a combina¬ 
tion of factors, but most importantly the ratio between the 
FFP’s initial ejection velocity and the escape velocity at the 
initial location of the FFP. Since we set up the FFP plan¬ 
ets following a Plummer distribution, the initial cumulative 
number of planets N(r) as a function of distance r to the 
cluster centre is given by: 


/ 2 \ — 3/2 

iV(r)=^l + ^J 


( 8 ) 


where N is the total number of planets in the system. The 
number of FFPs that escape at early times can then be 
obtained by combining Eqs. [5] and [8] and integrating over 
the entire cluster. Escapers are removed when they reach a 
distance of 2 rt (Eq. [2j|. For high-velocity FFPs this escape 
time is roughly 2rt(M )/ ctffp « 1 — 10 Myr. 

It is possible to estimate the evolution of the FFPs be¬ 
yond the initial phase of escape analytically. The linear rela¬ 
tion between f bs and fh P in Figure [4] is remarkable, and can 
be explained as follows. The number of stellar members of 
the cluster decreases with time, and this trend can to first 
order be approximated as linear: 


N sh (t) = N sb (0) (l- -M 

\ ^diss / 


(9) 


(e.g., Baumgardt & Makino 2003 Heggie & Hut 2003 
Lamers et al. 1120*05 1. With this simple relation, the time at 


which half of the stars have escaped from the star cluster is 
t ha = f d iss/2. The planet-to star ratio 77(f) = 7V p (f)/A r sb (t) 
decreases rapidly in a certain interval 0 < f < fi as high- 
velocity planets escape. Beyond this time, a certain fraction 
x = Ap(fi)/A7p(0) of the FFPs remains. The expression for 
77(f) for f > fi is thus 


77(f) 


TVp(fl) / fdiss ~ t \ 
7V sb (fl) V fdiss ) 


( 10 ) 


The number of FFPs in the star cluster as a function of time 
can then be computed using Eqs. 0 and [lO] 


JVp(f) = 77(f)7V sb (f) 


^p(*l)W sb (0) , 

f ‘HP 

Vi 0 

JVsbdi) 

^ tdlBI 

3-tl J tdlSB J 


= iV p(*i)(frfl) 


(ID 


The time f bp at which half of the planets have escaped from 
the star cluster can then be obtained by solving Eq. [Tl] for 
7V p (f)/7Vp(0) = 1/2: 



( 12 ) 


where x = N p (ti)/N p (0) indicates the fraction of FFPs that 
remain bound to the star cluster during the time interval 
[0, ti]. The above equation is valid for N p (ti)/N p (0) > 1/2, 
as larger initial escape fractions result in fh p < fi. In the 
limit where all FFPs are initially virialized, so that x = 1 
and fi = 0, Eq. |12| reduces to 

f hp = (2 - \/2)ths « 0.59fhs (13) 

In other words, if the FFP population in a star cluster ini¬ 
tially has a velocity dispersion similar to the velocity disper¬ 
sion of the star cluster, the planets typically escape at 60% 
of the remaining lifetime of the star cluster. This result is 
independent of the star cluster parameters, and provides a 
good estimate of the fate of a planetary population, under 
the condition that Eq. [9] is a reasonable approximation. 

From Figure [4] we can estimate the fraction of planets 
(1 — x) that escapes during the initial time 0 < f < fi by 
solving Eq. [l2]for x: 
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^ _ 1 f 2ihs — ti \ ^ 1 / 2ths \ 

2 \ 2 ths — Fp ) 2 2tb s — Fp ) 

This estimate is only valid when all FFPs are generated at 
t = 0, which is generally not the case in realistic star clusters. 
However Eqs. |12| and |14| but it can also be used to statisti¬ 
cally estimate long a FFP typically remains a member of a 
star cluster. 

The results in Figure [4] are well described by Eq. |12| 
The relation between tb p and ths is indeed linear, although 
slightly steeper in the regime where the star clusters can 
retain the majority of the FFPs for a substantial amount of 
time. The relation breaks down when many FFPs escape at 
early times. The time th p depends strongly on the ejection 
velocity distribution of the FFPs. A comparison between the 
panels in Figure [4] also shows that the FFPs escape earlier 
from clusters with a smaller velocity dispersion (i.e., a larger 
rinn). For / e j = P (1M 0 ) the ejection velocity of the planets 
is largest, and therefore thp is small; many planets escape 
within 10 — 20 Myr. 

3.1.3 Escape of the planet population 

Due to the preferred ejection of FFPs from the system, the 
planet-to-star ratio 77 (t) = N p (t)/N s b(t) rapidly decreases 
over time. The evolution of lZ(t) is roughly linear for most 
of the time, and we therefore quantify this evolution with a 
linear fit: 

^) = SbW~^ /(0)[1 + a/tl ’ (15) 

where t is the time in units of Myr and a in units of Myr -1 . 
In the ideal case where this dependence is truly linear, the 
quantity 77/(0) is the initial planet-to-star-ratio, and the 
quantity a/ ss dlZ(t)/dt indicates the relative escape rate of 
planets in units of Myr -1 . Although this derivative varies 
slightly with time, it is to first order constant. Under this 
assumption, the timescale at which all planets are removed 
from the star cluster is tdiss ~ —aj 1 . Apart from statis¬ 
tical fluctuations, a/ is independent of the initial number 
of FFPs, N p . Since the number of planets decreases more 
quickly than the number of stars, a/ is always negative. The 
quantity a/ can be used to estimate the number of planets 
in the star cluster at any time. 

In the case of a FFP population that is initially close to 
virial equilibrium with the stellar population, which is the 
case in Figure [2] the fitted value 77/(0) is close to (but not 
equal to) the initial planet-to-star ratio. When planets are 
ejected from their host systems with high velocities, a large 
number escape immediately, and those in the low-velocity 
tail of the distribution follow Eq. [15] Due to early escape, 
the fitted value 77/(0) is smaller than the initial planet-to- 
star ratio 77(0). This difference is thus an indicator of the 
initial fraction of high-velocity FFPs. 

The fitted parameters 77/(0) and a?/ for the different 
models are plotted in Figure[3] The top panel shows the val¬ 
ues for the models with rh m = 0.38 pc and f e j = P (1 Mq). 
All models have —1.4 x 10 3 < a/ < —0.5 x 10 -3 , which 
indicates that all planets are ejected from the clusters on 
timescales of roughly 0.7 — 2 Gyr. For models with identical 
values of N the scatter in 77/(0) is small and concentrates 
around the initial planet-to-star ratio 77(0). The time scale 


tdiss at which all planets escape from the cluster increases 
with increasing N s b, which corresponds to the dependence 
of the relaxation time thm on 7V a b (Eq. |6|. In addition, the 
results are (apart from statistical differences) independent 
of the number of FFPs, N p . 

Figure [3 ]d shows a comparison for clusters with identical 
Np and f e j, but with a varying number of stars 7V a b (differ¬ 
ent symbols) and half-mass radii rhm (different colors). Ini¬ 
tially, 77(0) = Np/N a h varies between 0.33 and 2. The ejec¬ 
tion timescale for the FFPs ranges from tdiss ~ 700 Myr for 
N s b = 500 to tdiss = 2 Gyr for N s b = 3000. Larger values of 
r v i r result in a smaller fitted value 77/(0) as compared to the 
initial conditions, as more FFP speeds are initially above the 
escape velocity (Eq. [5| and therefore escape at early times. 
For the clusters with the largest number of particles (where 
most FFPs are initially below the escape velocity), a/ in¬ 
creases with increasing JV S b, indicating that clusters with a 
smaller radius eject their planets faster, due to their shorter 
relaxation time. 

Figure [3J: shows the dependence of 77/(0) for star clus¬ 
ters with identical rh m and N p , but with different numbers 
of stars N s b = N(1 + B) and ejection velocity distributions 
f e j. After this initial phase of escape of high-velocity FFPs, 
the ratio 77(t) decreases roughly linearly with time. As this 
initial phase is short, and as a is independent of N p beyond 
this time, the corresponding timescale for the escape of the 
entire planetary population is still well-approximated with 
tdi aa = — aj , irrespective of the choice for f e j. 

The time at which the planet-to-star ratio 77(f) drops 
to half the value it had after the initial rapid escape phase, 
is tdiss/2. We quantify this half-life time for the stellar pop¬ 
ulation with th a , and for the FFP population with th p . The 
relation between ths and th P for the different star clusters 
is shown in Figure [4] From top to bottom, the three panels 
show the results for models with rhm = 0.38 pc, 0.77 pc, 
and 1.54 pc, respectively. The initial number of planets in 
each cluster is N p = 1000. Since planets escape faster due to 
their initially supervirial state and/or two-body relaxation, 
all star clusters have th s > th P , and both quantities are inde¬ 
pendent of Np. Note that the choice for N p does not affect 
the results in the figure, as the planets effectively behave as 
test particles. 

3.2 Close encounters 

During the simulations we record all close (< 1000 AU) en¬ 
counters between the members of the star cluster. Two-body 
encounters can be of the type star-star (S-S), star-planet 
(S-P), and planet-planet (P-P). The effect of gravitational 
focusing for S-S encounters is largest, while it is almost neg¬ 
ligible for P-P. Three-body encounters occur regularly and 
almost exclusively involve stellar binary systems (SS-S and 
SS-P). Three-body encounters that result in the formation 
of a dynamical binary system almost always result in the 
ejection of the body with the lowest mass. The formation of 
an SS binary is therefore much more common than a binary 
of type SP (a planet orbiting a star), while those of type PP 
(binary planets) only rarely occur. Moreover, FFPs tend to 
escape earlier than stars. As a result, close encounters of 
type SP-S are more common than those of type SP-P, even 
when initially N s < N p . Three-body encounters involving 
more than one planet in a bound system (SP-P and PP- 
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S) are rare, since this first requires the formation of these 
systems by dynamical capture. As initially all planets in 
the star clusters are free-floatintQ the encounter frequen¬ 
cies for the interactions SS-SP, SS-PP, SP-PP and PP-PP 
are negligible, primarily because capture of FFPs into SP- 
type systems is rare, while double planets (PP) almost never 
form. The only four-body encounters that occur frequently 
are those between two stellar binary systems (SS-SS), par¬ 
ticularly when primordial binaries are present. Despite these 
arguments, it should be noted that planetary systems may 
be common in real star clusters, and our results are therefore 
lower limits for the encounter frequencies of planet-hosting 
stars in open clusters. 

The encounter frequencies at a given time can be esti¬ 
mated using the abundances of the different types of sys¬ 
tems. Under the assumption that all particles have a sim¬ 
ilar distribution in phase-space, and ignoring the effects of 
gravitational focusing, the number of encounters N of the 
different types can be estimated through the following pro¬ 
portionalities: 


A f s -s 

«w s 2 

= (1 -BfN 2 

Ws-p 

oc N S N P 

= (1 - B)NN P 

A fp~p 

oc Alp 


Msp-s 

oc N sb N s N p 

= {l-B 2 )N 2 N p 

A fss-s 

oc N s N h 

= (1 - B)BN 2 

Nss-p 

oc N b N p 

= BN N p 

A fss-ss 

oc JVb 

= B 2 N 2 


In addition, all encounter rates depend on the choice of 
the encounter criterion Sr < p m and the size of the star 
cluster r v i r through the approximate proportionality A f oc 
{jpm / r vir (t) ) 3 at time t. These estimates are valid for short 
periods of time where the populations do not change sub¬ 
stantially though the formation or destruction of gravita¬ 
tionally bound two-body systems or dynamical mass segre¬ 
gation. 

3.2.1 Encounter rates for dynamical populations 

The frequency of occurrence of the various types of encoun¬ 
ters for a subset of the modelled star clusters with different 
initial conditions are listed in Table [2] The quantities repre¬ 
sent averaged results for an ensemble of ten realisations of 
each model. Note these are cumulative values over the entire 
lifetime of the star clusters. The relations in Eq.|16| therefore 
only provide first-order estimates. Model 3 is the reference 
model that is used for comparison with models that have 
different initial conditions. 

Models 1 — 3 represent the results for star clusters with 
different N and identical initial planet-to-star ratios, and 
shows that the number of encounters A f grows faster than 
the number of particles in the star cluster. This is because 
these three clusters have the same r\, m and therefore a stel¬ 
lar density and a total lifetime that increase with N. For the 
high-density star clusters the ratio Mp-p/Ns-s is substan¬ 
tially larger than for low-density star clusters. The reason 

1 Close encounters involving star-planet system (SP) may be 
common in realistic star clusters as a large fraction of stars is 
likely to host one or more planets. In our simulations, however, 
all stars are initialised without planetary companions. 


for this is that the ejection velocity for the planets is iden¬ 
tical for all models, and therefore planets in the low-mass 
star clusters are supervirial and escape at earlier times. Al¬ 
though the initial binary fraction is zero for these three mod¬ 
els, several binary systems (SS and SP) form and have close 
encounters with the other members of the star clusters. 

Models 3 — 5 show the results for models with an identi¬ 
cal N = N s + Nb and N p , but with different binary fractions 
B. The total number of individual stars in the clusters in¬ 
creases as N s b = N s + 2Nb = N{l + B), and the total cluster 
mass follows the same proportionality. The number of S-S 
and S-P encounters initially decreases with increasing B, but 
due to higher mass density, the destruction of binary systems 
and the longer lifetime of star clusters with a higher B, the 
cumulative number of these types of encounters increases 
with increasing B. For similar reasons, the P-P encounter 
rate, which is initially independent of B, increases with B 
due to dependence of the dissolution timescale on B. The 
total number of SS-S, SS-P, and SS-SS encounters increase 
strongly with B, although the occasional formation and de¬ 
struction of binary systems changes the proportionality of 
Eq. [16] mildly. 

The dependence of the the encounter frequencies on the 
number of planets N p can by studied by comparing mod¬ 
els 3, 6, and 7. As expected, the encounters involving only 
stars (S-S, SS-S, SS-SS) are independent of N p . The num¬ 
ber of encounters between stars and planets (S-P, SS-P) are 
proportional to N p , while for the planet-planet encounters 
(P-P) the proportionality is N p . The number of SP-S en¬ 
counters is well described by Eq. [16] with proportionality 
constant 7.5 x 10 -10 . 

Models 3, 8 and 9 represent identical star clusters apart 
from their initial size r v i r . Since the initial virial radius of 
model 8 is twice as large as as that of model 3, we expect all 
encounters involving only stars (S-S, SS-S, SS-SS) to occur 
roughly eight times less frequently. The number of encoun¬ 
ters between FFPs occur even less frequently because the 
ratio between the planet ejection velocities and the stellar 
velocity dispersion becomes larger as r v i r increases, resulting 
in the initial escape of a larger number of FFPs. 

Models 3 and 10 — 13 represent five star clusters with 
(statistically) identical stellar populations, but with differ¬ 
ent planet ejection velocity distributions / e j. Model 10 rep¬ 
resents the special case where all planets are assumed to be 
ejected from their host star with zero velocity, i.e., the FFP 
population in this model is initially in virial equilibrium. As 
expected, the number of encounters of types S-S, SS-S, SS- 
SS are statistically identical. All encounters involving FFPs 
depend on the ability of the star cluster to retain the escap¬ 
ing planetary population. The number of encounters involv¬ 
ing FFPs (S-P, P-P, SS-P) increases when the clusters can 
retain the FFPs for longer times, such as when the typical 
ejection velocity (Figure [T]) is lower. The number of FFPs 
that are captured into SP-systems is largest when the ejec¬ 
tion velocities are smallest. Models in which the FFPs have 
lower initial velocities therefore result in a larger number of 
encounters of type SP-S. 

3.2.2 The encounter number distribution 

In our simulations we also record the total number of close 
encounters experienced by each individual object. This num- 
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Table 2. The frequency of the various types of encounters during the lifetime of a subset of the star clusters. All values and their errors 
represent the average and the standard deviation for an ensemble of ten realisations. Several types of encounters are very rare (such as 
SS-SP), and are therefore not listed here. 


ID N 

B 

% 

N p 

r v ir 

pc 

M s 

Mq 

/ej 

S-S 

(1) 

S-P 

(2) 

P-P 

(3) 

SS-S 

(4) 

SP-S 

(5) 

SS-P 

(6) 

SS-SS 

(10) 

1 

500 

0 

1000 

0.5 

1.0 

p 

872 ±119 

357 ± 183 

77 ± 83 

93 ± 15 

1± 1 

7 ± 5 

2 ± 2 

2 

1000 

0 

2000 

0.5 

1.0 

p 

4374 ± 461 

2590± 1139 

656 ± 648 

254 ± 31 

2 ± 2 

17 ±8 

7 ± 3 

3 

2000 

0 

4000 

0.5 

1.0 

p 

22880 ± 2105 

20316 ± 6745 

6585 ± 5374 

707 ± 87 

7 ± 5 

57 ± 17 

13 ±4 

4 

2000 

20 4000 

0.5 

1.0 

p 

26265 ± 2044 

22959 ± 7741 

7700 ± 5901 

7599 ± 796 

29 ± 13 

2848 ± 956 

587 ± 100 

5 

2000 

50 4000 

0.5 

1.0 

p 

33814 ± 3462 

28254 ± 7824 

9523 ± 6539 

18817 ± 1555 

90 ± 23 

7472 ± 2144 

2715 ±277 

6 

2000 

0 

1000 

0.5 

1.0 

p 

21931 ± 1350 

4723 ± 1450 

377 ± 287 

726 ± 72 

1± 1 

16 ±7 

13 ±4 

7 

2000 

0 

2000 

0.5 

1.0 

p 

22458 ± 2016 

9796 ± 3265 

1602 ± 1296 

736 ± 86 

3 ± 2 

33 ± 16 

11 ±6 

8 

2000 

0 

4000 

1.0 

1.0 

p 

8190 ± 745 

4548 ± 2073 

1169 ± 1232 

457 ± 56 

1± 1 

28 ±9 

14 ±6 

9 

2000 

0 

4000 

2.0 

1.0 

p 

3286 ± 385 

1051 ± 593 

190 ± 241 

255 ± 34 

1± 1 

9 ± 9 

7 ± 3 

10 

2000 

0 

4000 

0.5 

- 

z 

21916 ± 1826 

42707 ± 2314 

27029 ± 1348 

695 ± 78 

15 ±5 

139 ± 24 

11 ±4 

11 

2000 

0 

4000 

0.5 

1.0 

D 

21875 ± 1670 

36927 ± 7276 

20778 ± 5758 

694 ± 77 

10 ±7 

119 ±29 

13 ±4 

12 

2000 

0 

4000 

0.5 

0.6 

D 

22070 ± 2225 

36351 ± 3993 

19765 ± 3203 

726 ± 51 

12 ±6 

126 ± 17 

12 ±4 

13 

2000 

0 

4000 

0.5 

0.6 

P 

19951 ± 7343 

25574 ± 9194 

10612 ± 3806 

627 ±227 

8 ± 5 

98 ±36 

9 ± 5 






Figure 5. The distribution of the total number of close (< 1000 AU) encounters per star for model 3 (top) and model 5 (bottom) 
during the entire lifetime of the star clusters. Both star cluster models initially contain N = N s + N b = 2000 (single or binary) stellar 
members with a binary fraction B = 0% (model 3) or B = 50% (model 5). Both have an initial planet-to-star ratio 7Z = 2 and an initial 
half-mass radius rhm — 0.38 pc. Encounter distributions are shown for encounters experienced by single stars (left) and by FFPs (right). 
The different colours indicate different types of encounters. The horizontal axis indicates how many encounters are experienced by a 
cluster member, and the vertical axis represents the number of cluster members that experience that number of encounters. Each of the 
histograms represents the average of an ensemble of ten realisations. 
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Figure 6. The number of FFPs that never approach a star or 
a binary system within 1000 AU during the entire lifetime of 
the star cluster. The horizontal axis represents the different mod¬ 
els, arranged by the initial number density in the star cluster. 
Each model initially contains N p = 1000 free-floating planets. 
The adopted initial planet ejection velocity distributions include 
prompt ejection (P), delayed ejection ( D ), and zero ejection ve¬ 
locity (Z). and are indicated with the different symbols (see Fig¬ 
ure [IJ. FFPs our lowest-density models (N sb = N s + N b = 500; 
B = 0%; r v ; r = 2 pc) almost never encounter a star or bi¬ 
nary system, while many FFPs in the highest-density clusters 
(N ab = 2000; B = 50%; r v ; r = 0.5 pc) do approach one or more 
massive cluster members within 1000 AU. In the latter clusters, 
the FFPs almost all FFPs experience a close encounter for the 
model in which the planet ejection velocity is zero (/ e j = Z). 


ber depends strongly on the location of the object in the star 
cluster and on the time it spends in the cluster before escap¬ 
ing. The (combined) total number of encounters experienced 
by all the objects of type S, P, SS, and SP is: 


A/"s,tot = 2A/"s-s + A fs-p + Wsp-s + A fss-s 

A/p,tot = A fs-p + 2Afp~p + J\fss-p 
A/"ss,tot = A fss-s + Wss-p + 2JVss-ss 
A/sp,tot = A fsp-s 


(17) 


where we have ignored the encounter types that rarely oc¬ 
cur, such as A fsp-p, Nsp-sp and A fss-sp (see also be¬ 
low). Note that during an encounter between two objects of 
the same type, both objects experience an encounter. Us¬ 
ing Eqs. |16[ and under the assumption that gravitational 
focusing is negligible (which may only be a reasonable as¬ 
sumption for encounters involving low-mass bodies), these 
can be expressed as 

A/"s,tot ocAi 2 (l-e)[2-e-|-A/ p (Ai- 1 + l-|-Z3)] 

A/p,tot <xN p (2N p + N) nsl 

A/ss,tot oc N 2 B( 1 + Np/N) 1 j 

A/sp,tot oc (1 - B 2 )N 2 N p 

at an instantaneous time f, under the assumption that the 
different populations have the same distribution in phase- 
space. In the remaining part of this section we discuss the 
encounters experienced during the entire lifetime of the star 
clusters, and since all variables in Eq. [18] change with time, 
these expressions only provide first-order approximations. 

The distributions of the number of close encounters per 
star and per planet, split up into the different components 


in Eq . |17| are shown in Figure [5j for models 3 and 5 (cf. 
Table |2h . The left-hand and right-hand panels show the en¬ 
counters experienced by the stars and by the planets, re¬ 
spectively. Note that in this figure we show the number of 
encounters experienced by each body (instead of the number 
of two-body encounters occurring); each encounter corre¬ 
sponds to two individual bodies experiencing an encounter. 
The area under each histogram corresponds to the values in 
Table [2] when the encountering cluster members are of dif¬ 
ferent type (e.g., S-P), and are double the values in Table[2] 
when the encountering members are of the same type (e.g., 
P-P). 


The top-left panel of Figure [5] shows the distribution 
of the number of encounters experienced by single stars in 
model 3. Based on our initial conditions, we expect to first 
order that the total number of encounters of each type is 
2A fs-s oc 2N a , JVs-p oc N p , and JVss-s ~ 0, for encoun¬ 
ters with other single stars, with FFPs, and with binary 
systems, respectively. Note that for our choice of the ini¬ 
tial conditions, we expect the distributions for the star-star 
and star-planet encounters to be roughly the same. However, 
due to the large fraction of FFP escapers at the early stage, 
the evolution of the star cluster, the formation of new dy¬ 
namical binaries, and gravitational focusing, the final results 
differ. Approximately 70 stars (3.5% of the total number 
of single stars) experience only one encounter with another 
star, and a slightly higher number of single stars experiences 
only one encounter with a FFP. Several single stars, on the 
other hand, experience over a hundred encounters with other 
stars, while none of the single stars encounters more than a 
hundred FFPs. Close encounters of single stars with binary 
systems are substantially less frequent, as the initial binary 
fraction is zero and only few binary systems are formed at 
later times. Approximately 3.5% of the single stars experi¬ 
ence close encounters with one or more dynamically formed 
binary systems, and most of these occur after the clusters 
have over half of their stars. 

For model 5 (bottom-left panel in Figure[5|, the number 
of single stars with a close encounter with one other single 
star is about 100 (10%). Although one would expect this 
number to be lower due to the smaller number of single stars 
in this model, this value can be explained as a result of the 
higher mass density in the star cluster, the disruption of wide 


binaries, particularly in the cluster centre (see, e.g., de Grijs 
et al.|2013), and the longer cluster lifetime, as compared to 


model 3. Approximately 130 single stars (13%) experience 
a close encounter with a FFP. The number of single stars 
experiencing a close encounter with a binary system is now 
substantial, due to the initial condition B = 50%. Single 
stars having many encounters with binary systems occur 
less frequently compared to single-single encounters, which 
is partially due to the disruption of the widest binaries over 
time, and partially due to the fact that binary systems tend 
to increase the velocities of single stars after an encounter, 
such that they are more likely to escape the star cluster. 

The distribution of close encounters experienced by the 
planets in models 3 and 5 are shown in the right-hand panels 
of Figure [5] For model 3, approximately 400 (10%) of the 
FFPs experience encounters with only one FFP, and a lower 
fraction (5%) of the FFPs have only one encounter with a 
star. All FFPs in the cluster experience fewer than 60 en¬ 
counters with another star, and fewer than 30 encounters 
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with another FFP. FFPs experience on average more en¬ 
counters with single stars and with FFPs in model 5 (w.r.t. 
model 3), since the number density for this model is larger. 
Again the encounter distributions for encounters with single 
stars and with FFPs are different. The difference in S-P and 
P-P encounters between the two models is partially a re¬ 
sult of the higher mass concentration, and partially a result 
of the longer lifetime of model 5. Close encounters between 
stars and binaries do not occur frequently in model 3 be¬ 
cause of the zero initial binary fraction, but as dynamical 
binaries form, several hundreds of FFPs have a close en¬ 
counter with a binary system. Encounters of the type SS-P 
are more frequent in model 5, and the majority of the FFPs 
have one or more encounters with a binary system. Although 
the initial binary fraction in model 5 is 50%, fewer FFPs ex¬ 
perience multiple close encounters with binaries than with 
single stars. 

The encounter distributions for the different types of 
cluster members in all models have in common that few bod¬ 
ies experience many encounters, while many bodies experi¬ 
ence few (or no) encounters. The former are mostly bodies 
that remain part of the star clusters until it is nearly dis¬ 
solved, while the latter escape or migrate the the cluster 
outskirts at early times. This means that statistical aver¬ 
ages, such as the results in Eqs. [Tt] and [18] normalised by 
the number of bodies, should be avoided when estimating, 
for example, the typical effect of close encounters on existing 
planetary systems. 

The fraction of FFPs that never experience a close en¬ 
counter with a stellar-mass particle (single or binary) is 
shown in Figure [6] for the different initial conditions. Note 
that the initial number density along the horizontal axis is 
proportional to the initial mass density, as we adopt the 
same IMF in all models. Thus this result also represent 
the correlation between FFPs with no encounters and clus¬ 
ter initial mass density. In the lowest-density star clusters 
none of the FFPs experience a close encounter, while in the 
higher-density star clusters the fraction depends strongly 
on the initial velocity distribution of the FFPs. The FFP 
populations that are initially in virial equilibrium (the filled 
circles in Figure [6] provide lower limits, and for these pop¬ 
ulations the fraction of non-encountering FFPs ranges from 
0% to 100%, primarily depending on the star cluster density. 


3.2.3 Encounter times 

The cumulative number of encounters Mx-y{t ) between 
members of two types of populations X and Y should follow 

Ux-y oc f Nx(t)NY(t)dt (19) 

Jo 

where Nx and Ny are the number densities of these popu¬ 
lations, respectively. When making the assumption that the 
proportionality remains constant (i.e., when we ignore mass 
loss, mass segregation, core collapse, and cluster expansion), 
the time th at which half of the encounters have occurred 
can be obtained by solving A fx-Y^th) = A/"x-y(fdiss)/2 for 
th- Approximate values for th for the various types of en¬ 
counters can be derived by substituting Eqs. [9] and/or |11| 
into Eq. [19] giving 





Figure 7. The times at which 10% (top), 50% (middle), and 
90% ( bottom ) of the encounters have occurred, for all models. The 
horizontal axis represents the initial number density in the star 
cluster, from the lowest-density models (A r s |, s= IV S + IV), = 500; 
B = 0%; r vir = 2 pc) on the left to the highest-density clusters 
(N s b = 2000; B = 50%; r v ; r = 0.5 pc) on the right. The vertical 
axis represents the time at which a certain fraction of encoun¬ 
ters have occurred (cf. Eq. |20| ). Different colours represent differ¬ 
ent encounter types between stars (S), FFPs (P), binary systems 
(SS) and captured planetary systems (SP). The adopted initial 
planet ejection velocity distributions include prompt ejection (P), 
delayed ejection ( D ), and zero ejection velocity (Z), and are in¬ 
dicated with the different markers (see also Figure [TJ . The data 
fl.rp Hp.rivprl from an p.nspmhlp of t.p.n rfiA.lisat,inns. 
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th,s-s « (1 — 2 1 ' /3 )tdiss ~ 0.21idiss 

th,s-p « (1 - 2- 1/4 )fdis S ~ 0.16i d iss (20) 

th,P-P ~ (1 — 2 -1 / 5 )tdiss «0.13tdiss 

where we have for simplicity assumed that the FFP pop¬ 
ulation is initially virialized, that B = 0, and that Eq. [9] 
holds . The expressions can be generalized for the cases with 
ATp(ti) < JVp(O) and t\ > 0 using Eq. [lOj and will result in 
smaller values of M{th) for the encounters involving planets. 
Note that Eq. [20]is independent of the number of stars and 
the number of planets, if the assumptions stated above hold. 

The timescales at which 10%, 50% and 90% of the en¬ 
counters of the specific types occur, are shown in Figure [7] 
for all modelled star clusters. Note that we model each star 
cluster until it has completely dissolved. The times at which 
these events occur vary with the initial properties of the star 
clusters and their FFP populations, and also with the types 
of encounters. In general, however, 10% of the encounters 
occur in the first 3 Myr, 50% occur within 30 Myr, and 90% 
of the encounters occur within 100 Myr. 

The middle panel shows the times at which half of the 
encounters have occurred, and its values can be approxi¬ 
mated with Eq. |20| Although tdiss depends on the star clus¬ 
ter properties (particularly on r v j r , N s b and B), the ratios 
th,s-s/th,s-p ~ 1-31 and t h , S -s/th,P-p « 1.62 are inde¬ 
pendent of these, which explains why the differences in th in 
Figure [7] are roughly constant in logarithmic units. The real 
differences in the figure are slightly larger than in Eq. [20] 
because several of the assumptions are violated, as a result 
of expansion and early loss of FFPs in the systems. This can 
be observed in the top panel, which shows that the 10% en¬ 
counter time of the P-P interactions is substantially smaller 
than the simple approximation suggests for the models with 
high planet ejection velocity. For these cases, a better ex¬ 
pression may be derived by not adopting N p {ti) = N p {0) 
and t\ = 0 in Eq. m as done above. However, even though 
a large number of planets may escape at early times, the ma¬ 
jority of the P-P encounters occur between the FFPs that 
remain in the star clusters for a long time, and therefore the 
results in the middle and bottom panels depend only mildly 
on / e j. The encounters involving massive components (such 
as SS-S and SS-SS) generally occur at later times than those 
with low-mass components, as the former remain part of the 
star cluster for longer times. Figure[7]also demonstrates that 
encounters occur at later times for systems with increasing 
N s b, B, or r v i r , since these systems have larger relaxation 
times (Eq. [6| and generally have longer dissolution times. 


3.3 Encounter velocities and periastron distances 

3.3.1 Encounter velocity distributions 


The trajectories of two approaching unbound bodies in iso¬ 
lation are hyperbolic, and have an associated velocity-at- 
inhnity Voo- The velocity-at-infinity distribution, f(voo), of 
an ensemble of two-body encounters in a star cluster is ap¬ 
proximately a Maxwell-Boltzmann distribution: 

-^r') 47 tv;L (21) 



where v-m is the most probable value of Voo . The correspond- 
ing mean value is (voo) = 2n~ 1 ' 2 v m . Eq. 21 is generally a 


good approximation and is independent of mass of the bod¬ 
ies involved in the encounter, provided that the cluster is not 
mass segregated and that the star cluster’s velocity distri¬ 
bution is isotropic. For an encounter with impact parameter 
b, the velocity at periastron can be expressed as 

with e = \/ 1+ (§z) ’ (22) 



where e is the hyperbolic eccentricity. The distribution of 
impact parameters follows f(b) oc b. When the combined 
mass of the bodies is small, such as for P-P encounters, 
the equation reduces to v p = Voo- For a close, low-Uoo en¬ 
counter between massive bodies, the expression reduces to 
v p « 2GM/bvoo. We approximate the distribution of the 
velocities at periastron, f{v p ) with a modified Maxwell- 
Boltzmann distribution that has a velocity offset v 0 and a 
dispersion parameter v m : 

^ An(v p - Vo) 2 , (23) 


f(v p ) oc (irv^) 3/2 exp 


(v p - v 0 


where v 0 = v p — Voo > 0 is the velocity increase during the 
encounter and v m is the most probable periastron velocity. 

The histograms in Figure [8] show the distributions of 
Voo {top) and v p {bottom) for the various types of encoun¬ 
ters that occur in star cluster models 3 {left-hand panels) 
and 5 {right-hand panels), and the corresponding best fits 
to Eqs. |21| and |23| The curves represent the distributions in¬ 
tegrated over all locations in the star clusters, for their entire 
evolution. Since the encounter velocities depend both on the 
location within the star cluster and on time, Eqs. |21| and |23| 
are approximations. Nevertheless, since most encounters oc¬ 
cur within the first ~ 100 — 200 Myr, the approximations 
are reasonable. 

The strongest deviations from the Maxwellian distribu¬ 
tion are seen for the largest values of Voo and v p , indicating 
interactions in the centre of the star cluster, and for the 
smallest values of Voo and v p , indicating encounters in the 
cluster outskirts and during the latest phases of the evolu¬ 
tion when the stellar density is low. Since planetary systems 
have v p « Voo, the corresponding curves for the P-P en¬ 
counters are very similar in all four panels, and well-fitted 
by Eq. 21 All encounters involving planets (P-P, S-P, SS-P, 
and P-P) are qualitatively very similar, both for model 3 
and model 5. As a result of gravitational focusing, the typi¬ 
cal velocity at periastron increases as the total mass of the 
encounters increases. 

With respect to the Maxwellian fit, excess S-S inter¬ 
actions with small Voo occur. These interactions are most 
frequent in the low-density regions in the outskirts of the 
clusters and also at late stages of evolution, when the ef¬ 
fect of gravitational focusing on low-velocity stars is more 
prominent. The deviation from a Maxwellian distribution is 
much stronger for the SS-S interactions. Since model 3 does 
not have primordial binaries, all binaries involved in the SS- 
S encounters are formed through capture. As these binaries 
are generally formed in the cluster centre, and as their com¬ 
ponents are generally of high mass, Eq.[2l]does not provide a 
good approximation for the distribution. Since B = 50% for 
model 5, the effect of gravitational focusing and longevity of 
the binary population is now much better visible. The SS- 
S and SS-SS show deviations from the Maxwell-Bolzmann 
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(a) N = 2000, B = 0% TZ = 2, r hm = 0.38 pc 



(b) N = 2000, B = 50% K = 2, r hm = 0.38 pc 



(c) N = 2000, B = 0%, TZ = 2, r hm = 0.38 pc 



Figure 8. The velocity-at-infinity (v rK: ) distributions of the close encounters, for model 3 ( top-left ) and model 5 ( top-right ), and the 
periastron velocity (v p ) distributions of the close encounters for model 3 ( bottom-left ) and model 5 ( bottom-right ). The histograms 
represent the simulation results, and the curves are best fits to Eq. |21| for Voo and to Eq. |23| for v p . The data represent the combined 
results for ten realisations. 


distributions, because of gravitational focusing and because 
they tend to remain longer in the star cluster. Note that 
the dynamically-formed binary systems in model 3 have on 
average substantially larger v p than the primordial binaries 
in model 5, although they are fewer in number. As these 
dynamical binaries in model 3 are generally more massive 
than the primordial binaries in model 5, their encounters 
with other stars and FFPs are stronger, and therefore their 
periastron velocities are larger. 

There are about ten times more SP-S encounters in 
model 5 than in model 3, indicating that the presence of bi¬ 
nary stars in star clusters enhances the star-planet captures. 
The distribution f(v oo) for SP-S encounters is a scaled-down 
version of that of the S-S encounters, which may be expected 
since the planetary companions do not play a role in the en¬ 
counter trajectory or encounter rates. 

The best-fitting values of v rn and v„ in Eqs. |21| and |23| 
are plotted in Figure [9] for the different types of encounters, 
for all simulated models. Note, however, that we did not 
include the bad fits, which result either from distributions 
that are substantially different from Eq. |21| or from low- 
number statistics. 


Most encounters occur during the early phase of star 
cluster evolution when mass segregation is small, and dur¬ 
ing this period we expect v m to be directly proportional 
to the velocity dispersion of the star cluster (Eq. [I]). For 
star cluster members near the half-mass radius v m oc <r(r = 
rhm) oc (A a b/Vh m ) 1 ' /2 . This trend is clearly seen for most fits 
to f(v oo) plotted in Figure [9] and is independent of f e j, N p , 
and the type of encounters. Since the P-P encounters have 
v p « Voo , we obtain the same relation for f(v p ). For the en¬ 
counters between massive components the fitted values for 
v m are larger, and v m typically increases with increasing 
component masses. 

Since FFPs barely change their velocities during P-P 
encounters, we expect v 0 ~ 0. The scatter for the P-P inter¬ 
actions results from the fact that we have fitted the distri¬ 
butions f(v oo) and f(v p ) to all encounters integrated over 
the entire volume of the star cluster and its entire evolution. 
The value of v 0 depends only mildly on B, N sb , f e] and r hm . 
As the fitted value v m takes into account the dependency of 
f{v p ) on Voo, the value of v 0 is primarily determined by the 
mass of the encountering components. Since gravitational 
focusing is stronger for encounters between massive bodies, 
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increases with the combined 
mass of the members involved in the encounter. 


The fitting results for f(vp) are combined in Figure 10 


The value of v m is primarily determined by the velocity dis¬ 
persion of the star cluster, while v 0 is mostly determined by 
combined mass of the encountering bodies. These trends are 
seen in Figure fTO] models with identical velocity dispersion 
lie roughly on vertical lines, while the same types of encoun¬ 
ters (S-S, S-P, etc) lie roughly on a horizontal line. Apart 
from some scatter, the values of v 0 of most of the encounter 
types indeed lie roughly on a horizontal line in Figure [9j 
although there is a trend of v 0 increasing with v m (and 
therefore with increasing N a \,). This trend is strongest for 
the P-P encounters. The results in Figure [l0| can be used to 
estimate the periastron velocity distributions for other star 
cluster masses and for encounters between members with 
other masses, using Eq. |23| It can also be used to statis¬ 
tically estimate the probability of obtaining a certain pe¬ 
riastron velocity for a single free-floating planet in a star 
cluster. 


3.3.2 Encounter periastron distance distribution 


For hyperbolic orbits between two objects with masses Mi 
and M 2 , the relation between the impact parameter 6, the 
periastron distance p, and the velocity at infinity Voo is given 
by: 


b = p 



2 GM 

2 ’ 
pv So 


(24) 


where M = Mi + M 2 is the combined mass of the two ob¬ 
jects. For two approaching bodies in a virialized star cluster, 
the velocity at infinity Voo is typically the velocity disper¬ 
sion in the star cluster multiplied by a factor \/2. When both 
bodies have a negligible mass, such as in the case of a planet- 
planet encounter, the equation reduces to p = b, while for a 
close encounter between two large masses the effect of grav¬ 
itational focusing is important and the equation reduces to 
p = (bv oo) 2 (2 GM) 1 . As the distribution over impact pa¬ 
rameters scales as fb(b) oc 6, the corresponding distribution 
for p is f(p) oc p for distant encounters. The corresponding 
distribution for close encounters can be obtained with trans¬ 
formation of variables, which results in f(p) oc r^(2GM) _1 . 
In other words, f(p) is independent of p in this regime. 

Figure EH shows the periastron distance distributions 
f(p)dp for the different types of encounters for models 3 
and 5. The distributions are linear over most of the range in 
p, and we therefore describe the results using the functional 
form 


f(p)=H (ap + /3)= Af ( a(p - p m / 2) + p m J ) (25) 

to the distributions for the different types of encounters, 
where N represents the total number of encounters of a 
specific type, a is in units of AU -2 and /3 in units of 
AU” 1 . Normalization to the total number encounters gives 
ft = Pm 1 ~ apm/2., where p m = 1000 AU is the largest peri¬ 
astron distance that is taken into consideration. In the limit 
that the masses of the encountering bodies are negligible, 
f(p) = 2jVp/p^ n , while in the limit of very close encoun¬ 
ters between massive bodies, f(p) = /3 = Af /p m . The corre¬ 
sponding cumulative distribution is 





Figure 9. Properties of the velocity-at-infinity distributions 
f(v a o) and periastron velocity distributions f(v p ) for the differ¬ 
ent types of encounters for the modelled set of star clusters. The 
horizontal axis represents the initial stellar density, as in Figure[7] 
Top: the fitted parameter v m (Eq. |21| ) to the measured distribu¬ 
tion f(v 00 )• The two other panels show the best-fitting param¬ 
eters Vm ( middle ) and v 0 ( bottom ) in Eq. |23| for the measured 
distribution f(v p ). All models resemble model 3 (cf. Figure [sj. 
Different colours represent different encounter types. Data with 
fitting errors larger than 0.5 are omitted. Different markers in¬ 
dicate the initial planet ejection velocity distributions, including 
prompt ejection (P), delayed ejection ( D ), and zero ejection ve¬ 
locity (Z); cf. Figure [T|). Data are derived from an ensemble of 
ten model realisations. 
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Figure 10. Properties of the periastron velocity distributions 
f(v p ) for the different types of encounters for the modelled set of 
star clusters. The fitted parameters v m and v 0 to Eq. |23| are indi¬ 
cated along the horizontal and vertical axes, respectively, for the 
models with r^m = 0.38 pc, N p = 1000 and with ejection veloc¬ 
ity distributions corresponding to the prompt- IMq distribution 
in Figure^ The different markers indicate the different types of 
encounters between stars (S), FFPs (P), binary systems (SS), and 
captured planetary systems (SP). The values of N s b = N s + N b 
are indicated with the different colors. 


F( P )=AT[ Qp{p -Pm) + P_ 

2 Pm 


(26) 


which indicates the cumulative number of encounters smaller 
than p. Note that in this approximation the values of a (and 
/3) are independent of the number of encounters, and also in¬ 
dependent of N a b, N p , and B. For encounters between more 
massive bodies, however, a (and therefore also /3) depend 
on Upo, and the latter quantity may have a different distri¬ 
bution in space (due to mass segregation) and time (due 
to the preferred escape of low-mass members) than for P-P 
encounters. 

In realistic star clusters f(p) is a combination of weak, 
intermediate, and strong encounters, and for the combina¬ 
tion of all encounters of a certain type, integrated over the 
entire star cluster and over time. Nevertheless, f(p) is still 
reasonably well approximated by the expression in Eq. |25| 
Several of the distributions in Figure m show an overabun¬ 
dance of periastron approaches below 10 AU with respect to 
Eq. |25[ This overabundance is exclusively seen for encoun¬ 
ters involving single stars (S-S, S-P, and S-SS) in model 5, 
and is related to the evolution of primordial binary systems, 
their interactions with other star cluster members, and sub¬ 
sequent decay. 

Figure 


12 


shows the fitted values of a (Eq. 251 for all 


models. Since the vast majority of the P-P encounters are 
weak, the fitted value for a is close to the predicted value in 
the unfocused limit, a. = 2p“ 2 = 2 x 10 -6 AU~ 2 , although 
the scatter around this value is considerable, particularly for 
the low-density star clusters for which a large fraction of the 
FFPs escape immediately. Larger masses for the encounter¬ 
ing bodies generally result in smaller values of a, which is 
consistent with the increased gravitational focusing. In addi¬ 
tion, when the star cluster mass increases and/or its radius 




Figure 11. The periastron distance distributions f(p) of the close 
encounters for model 3 (top) and model 5 (bottom). Different 
colours represent different types of encounters between stars (S), 
FFPs (P), binary systems (SS) and captured planetary systems 
(SP). The solid histograms indicate the number of close encoun¬ 
ters in each bin, and the dashed lines are linear fits. Encounter 
types that occur less than 200 times (e.g., SS-SS) are omitted. 
The curves represent the combined results for ten realisations. 


decreases, the velocity dispersion (and hence Voo) becomes 
larger, such that gravitational focusing is less important and 
the measured values of a are larger. 

The typical time that two encountering members spend 
near periastron is At = 2 p/v p . For weak encounters this 
is At ss 2b/Voo and for strong encounters this is At « 
bVoo/GM. This timescale is important in studying the ef¬ 
fects of encounters when one of the two members, for exam¬ 
ple a single star, has a planetary system, as it describes, 
given a planet orbiting a star with a certain semi-major 
axis, whether the encounter is impulsive or adiabatic, and 
whether it is close or tidal (see the extensive study of 
Spurzem et al. [2009, for further analysis). 


4 CONCLUSIONS AND DISCUSSION 

Free-floating planets (FFPs) may be abundant in young star 
clusters, as close encounters between stars can destabilise 
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types of encounters are well fitted by distributions f(v p — v 0 ) 
similar to a Maxwell-Boltzmann distribution, with a velocity 
offset Vo- The most frequent velocity v m in this distribution 
is proportional to the velocity dispersion of the star cluster, 
and the offset v 0 is primarily determined by the combined 
mass of the encountering bodies of a certain type. 

(v) The distribution over periastron distances is linear, 
f(p) = ap + P for the encounters is linear over most of the 
values of p. In our case we have recorded all approaches 
within p m = 1000 AU. In that case, distant encounters have 
a = 2 x 10 -6 AU -2 and j3 = 0 AU -1 . In the case of close 
encounters, f(p) = /3 = Af/p-m, where A f is the total number 
of encounters and p m is the close encounter limit. Also when 
we combine the encounters for the entire life of the star 
clusters, the linear fit is good, except for very small distances 
(p < 10 AU), for which approaches are more common. 


Figure 12. The fitted proportionality parameter a (Eq. |2o[ l for 
the periastron distance distributions f(p) for the different star 
clusters. The horizontal axis represents the initial stellar density, 
as in Figure [7] The adopted initial planet ejection velocity dis¬ 
tributions are indicated with the different markers, and include 
prompt ejection (P), delayed ejection (D), and zero ejection ve¬ 
locity (Z); see Figurejl] The different types of encounters between 
stars (S), FFPs (P), binary systems (SS) and captured planetary 
systems (SP) are indicated with the colors. Results are shown for 
a combined ensemble of ten realisations of each model. 


planetary systems, which can result in direct or delayed ejec¬ 
tion of planets from their host planetary systems. In this ar¬ 
ticle we have presented a study of the dynamical evolution 
of FFP populations in various types of star clusters. We 
have carried out IV-body simulations to characterise how 
the FFP population evolves and we have studied the prop¬ 
erties of close (< 1000 AU) encounters between single stars, 
binary stars, and FFPs. Our results can be summarised as 
follows: 

(i) A certain fraction of the FFPs escape from their star 
clusters shortly after they are ejected from their host stars, 
while the remaining FFPs escape at much later times. Be¬ 
yond the time of early escape, the planet-to-star ratio de¬ 
creases linearly with time, and the times ths and th P at which 
respectively half of the stars and planets escape from the star 
clusters are related through the simple relation th P ~ 0.6ths- 
For individual planets with low ejection velocities this means 
that they are likely to escape from the star cluster at a 40% 
earlier time than their host star. 

(ii) Many FFPs ejected from their host star system at 
early times experience tens of close (< 1000 AU) encounters 
with other stars and FFPs in the star cluster before escaping 
from the cluster. The fraction of FFPs that leave the cluster 
without any close encounter increases with increasing initial 
velocity distribution and with decreasing stellar density. 

(iii) Typically half of the encounters of all types occur 
within 30 Myr, while only 10% of the encounters occur after 


Our study is intended to obtain a general overview of the dy¬ 
namical evolution of a FFP population in star clusters. We 
have made several assumptions that should be kept in mind 
when making realistic predictions for existing star clusters, 
or when interpreting observational data. First, a critical as¬ 
sumption we made is that the FFP space density is identi¬ 
cal to the stellar density. In reality, the ejection velocities of 
FFPs may depend in an intricate way on the position and 
velocity of the star at the moment when it is perturbed. This 
is especially the case when planets are ejected promptly after 
their host system experiences a close encounter. One would 
then expect planets to be preferably ejected when their host 
star is near the star cluster centre. In the case of delayed 
ejection, however, this may well be a good approximation, 
since the stars are well-mixed at the moment of planet ejec¬ 
tion. 

Second, in our initial model setup we have also made the 
crude assumption that all planets are initially free-floating. 
We also neglect any additional FFPs produced by the close 
encounters between stars in our model, as well as the possi¬ 
ble capture of FFPs by planetary systems. Although these 
is an unrealistic assumptions, particularly in the case of per¬ 
turbed planetary systems that may eject planets tens of mil¬ 


lions of years later (e.g., Malmberg et al. 2011 Hao et al. 
|2013[ ), it is a good approximation if most planet ejections 
occur at early times. Moreover, many of the analytical and 
computational results presented in Section [3] are easily scal¬ 
able to any value of N p . These can therefore also be used to 
evaluate the probabilities for close encounters and escape for 
individual planets (A p = 1). When convolved with an time- 
dependent FFP production rate, estimates for the dynamical 
behaviour of the entire population of FFP can subsequently 
be obtained. 

In a follow-up study we will analyse more realistic initial 
conditions, ideally by modelling the full A-body evolution of 
decaying multi-planet systems in young star clusters, poten- 


the first 100 Myr. The ratios of the timescales at which half 

tially using the AMUSE framework ( 

Portegies Zwart et al. 

of the number of star-star, star-planet, and planet-planet 

2013 Pclupessy et al. 2013 Cai et al. 

20151, Specifically, we 

encounters occur are t h ,s-s ■ th,s-p ■ th,p-p ~ 1 : 0.77 : 

will analyse the dynamical evolution of a FFP population 


0.62. 

(iv) The velocity-at-infinity distributions, f(v oo), of the 
encounters are well approximated with Maxwell-Boltzmann 
distributions. The periastron velocity distributions for all 


where the FFPs have a stronger preference to be generated 
in the cluster centre, and where they are ejected at appropri¬ 
ate times resulting from close encounters with other cluster 
members and the decay of perturbed multi-planet systems. 
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